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0^ ■ Abstract 



We apply a recent proposal for defining states and observables in 
quantum gravity to some simple models. First, we consider the toy 
model of a Klein-Gordon particle in an external potential in Minkowski 
spacetime and compare our proposal to the theory obtained by de- 
parametrizing with respect to a choice of time slicing prior to quan- 
^^ . tization. We show explicitly that the dynamics defined by the de- 

V-j ' parametrization approach depends upon the choice of time slicing. On 

the other hand, our proposal automatically yields a well defined dy- 
namics - manifestly independent of the choice of time slicing at in- 
termediate times - but there is a "memory effect" : After the poten- 
tial is turned off, the dynamics no longer returns to the standard, 
free particle dynamics. Next, we apply our proposal to the closed 
Robertson- Walker quantum cosmology with a homogeneous massless 
scalar field. We choose our time variable to be the size of the universe, 
so the only dynamical variable is the scalar field. It is shown that 
the resulting theory has the expected semi-classical behavior up to the 
classical turning point from expansion to contraction, i.e., given a clas- 
sical solution which expands for much longer than the Planck time, 
there is a quantum state whose dynamical evolution closely approx- 
imates this classical solution during the expanding phase. However, 



when the "time" takes a value larger than the classical maximum, the 
scalar field becomes "frozen" at the value it had when it entered the 
classically forbidden region. The Taub model, with and without a ho- 
mogeneous scalar field, also is studied, and similar behavior is found. 
In an Appendix, we derive the form of the Wheeler-DeWitt equation on 
minisuperspace for the Bianchi models by performing a proper quan- 
tum reduction of the momentum constraints; this equation differs from 
the usual form of the Wheeler-DeWitt equation, which is obtained by 
solving the momentum constraints classically, prior to quantization. 

PACS #: 04.60.-hn, 03.70.-Fk, 04.20.Cv 



1 Introduction 

The "problem of time" refers to the difficulties in defining a Hilbert space of states, observables, 
and a nontrivial notion of dynamics in a theory where the spacetime metric is itself a dynamical 
variable, so no background metrical or causal structure is present. In the canonical approach to 
quantum gravity, this difficulty manifests itself directly by the fact that ordinary time evolution 
corresponds to a diffeomorphism (i.e., a gauge transformation), and the Hamiltonian is thereby 
constrained to vanish [p. As a result, the "wavefunction of the universe" in quantum gravity 
does not depend on "time" . This necessitates a serious reconsideration of the structure and 
interpretation of the theory provided by a naive application of the usual canonical quantization 
prescription. 

Recently one us proposed Q a prescription for defining states, observables, and dynamics 
in quantum gravity by making use of the background metrical and causal structure which is 
available in superspace. The main purpose of this paper is to further elucidate this proposal by 
applying it to some simple models. In Sec. Q we briefly outline the proposal in the context of 
quantum cosmology. In Sec. y we study a toy model of massless scalar particle in an external 
potential in Minkowski spacetime - which shares some of the key features of the quantum- 
cosmological models - for the purpose of comparing the proposal with the "reduced theory" 
obtained by deparametrization with respect to a choice of time variable. We show explicitly 
in this model that the dynamics defined by the deparametrization approach depends upon the 
choice of time variable. In contrast, the dynamics of our proposal is well defined, but there is a 
"memory effect" , wherein the dynamics does not return to that of standard free particle theory 
after the potential is turned off. In Sec. Q we then turn to an analysis of some features of the 
proposal in the context of general homogeneous cosmologies with a scalar field. In Sec. |5| we 
apply the proposal to the closed Robertson- Walker cosmology with a homogeneous scalar field. 
We find that the semi-classical description is a good approximation right up to the classical 
turning point from expansion to contraction. We also find that after this stage, the scalar 
field - which is the only dynamical variable - is frozen, i.e., it does not evolve with "time", 
which is chosen to be the size of the universe. In Sec. y we repeat the analysis of Sec. y for 
the Taub universe with and without a homogeneous scalar field and find qualitatively similar 
results. Finally, in Appendix A, we derive the Wheeler-DeWitt equation for the Bianchi models 
by performing a quantum (rather than classical) reduction of the momentum constraints. 

2 The proposal 

First we briefly review the proposal of Ref. |Q] in the context of quantum cosmology. The 
dynamical variables in quantum cosmology are the homogeneous spatial metric hab and possible 
homogeneous matter fields together with the conjugate momenta of these variables. Let us 
denote the configuration variables (i.e., coordinates on minisuperspace) hy q {A = 1, . . . ,N) 
and their conjugate momenta by tta- Then, as will be discussed in detail in Sec. ^ (see also 
Appendix A), for the models under consideration here, the classical Hamiltonian constraint 
takes the form 

H = G^^7TATTB + V{q)=0. (2.1) 

Here the supermetric. Gab, has Lorentz signature, — -|-- • • -|-, and, in our models, minisuperspace 
is globally hyperbolic in this metric. Furthermore, Gab possesses a one-parameter group of 
timelike conformal isometries corresponding to scaling transformations of the spatial metric. 



I.e., 

.2a 



hah -^ e'^hab . (2.2) 

(This one-parameter group of conformal isometries extends to full superspace Q, and, thus, 
is not an artifact of the simplicity of our models.) For convenience, we shall redefine Gab, 
if necessary, by a conformal transformation so as to make these conformal isometries be true 
isometries. 

We take the quantum version of the Hamiltonian constraint - i.e., the Wheeler-DeWitt 
equation - to be a Klein-Gordon equation of the form 

-V^VA + e^ + F(g)]$ = 0. (2.3) 

Since Gab is naturally defined only up to conformal equivalence class, it would appear most 
natural to choose ^ so that ( |2.3D is conformally invariant Q (see, however, the discussion 
following Eq. ([A9|) of Appendix A) . 

To define a Hilbert space structure on an appropriate subspace of solutions, $, of this 
equation, we make use of the natural (real) symplectic product, Q, given by 

n{^u ^2) = / [$2Va$i - $iVa$2] dS^ , (2.4) 

where S is any Cauchy surface. Note that $7 is conformally invariant ||5[ provided that we 
scale the wavefunction <^ as <I> ^ A(q')^~^'^$ under Gab — > ^{q)'^Gab- With our choice of 
supermetric, minisuperspace is static as a spacetime under the isometries defined by Eq. (|2.2|) , 
so if V{q) were a-independent, then the solutions that are positive frequency with respect to a 
could be chosen to form the Hilbert space [y], with inner product 

{^l,<^2)KG = -i^i^l,<^2). (2.5) 

However, the potential V{q) depends on a and, in fact, there is no symmetry of the super- 
Hamiltonian, H, that can be used to define positive frequency solutions [H]. Nevertheless, we 



have y— >Oasa— >— 00. Thus, the Wheeler-DeWitt equation (2^) possesses an asymptotic 
symmetry at "early times". It was conjectured in Ref. [^ that - at least in a suitable class 
of homogeneous cosmologies - this asymptotic symmetry enables one to construct the Hilbert 
space, TC, of solutions that are positive frequency with respect to a in the limit a -^ —00, with 
inner product (2^). We shall discuss this issue further in Sec. ^ 



Assuming that the Hilbert space, Ti, is obtained as we described, we complete our construc- 



tion of a quantum theory corresponding to the classical super-Hamiltonian (2^) by specifying 
the self-adjoint operators on TC which represent "position and momentum" observables at a 
given "instant of time". We take the notion corresponding to an "instant of time" in this 
theory to be the specification of a Cauchy hypersurface, S, in minisuperspace. There exists a 
standard prescription for the construction of the desired position and momentum operators on 
the Hilbert space ^^(S) (see, e.g., appendix C of Ref. 0) but our inner product is not the L"^ 
inner product but the Klein-Gordon one. Nevertheless, any C^ solution to the Klein-Gordon 
equation which lies in our one-particle Hilbert space Ti. is uniquely determined by its value 



on S 1^, since, by Eqs. ( p^ ) and (2^), the difference between two solutions in TC having the 
same restriction to S must have vanishing Klein-Gordon norm. As in Ref. Q, we assume that 
the subspace, V, of C^ solutions in TC whose restriction to S lies in L^(S) is dense both as a 
subspace of TC and as a subspace of L^(E), and that if a sequence in V converges in both TC and 
L^(S), then its limit in TC is nonzero if and only if its limit in -L^(S) is nonzero. If we view the 



Klein-Gordon inner product ( p.5| ) as a quadratic form defined on a dense domain in L^($]), it 
then follows that one can write 

(^1, ^2)a'G = (5($i|s), S($2|e))l2 , (2.6) 

where B : i^(S) — > L^(S) is a positive self-adjoint operator. We then define the position and 
momentum observables at "time" S by 

{<^i\0\<^2)kg = (5(«5i|e),05($2|e))l2 , (2.7) 

where O is the standard representation of the observable on L^(S). Note that the operator B 
defined by Eq. ( |2.6| ) depends, in general, on the entire "history" of the minisuperspace from the 
"asymptotic past" (i.e., a — > — oo) to "time" E, so the definition of O is nonlocal in "time" as 
well as "space" . 

The time evolution of a state ^ € Ti is, of course, given simply by Eq. ( |2.3D , but the 
observables, O, are not represented in a simple manner on TC. It is useful, therefore, to derive 
the evolution equation for the corresponding states ^ = i?(<I>|s) in L^(S), where the observables, 
O, take a simple form. Let S^ be a one-parameter family of Cauchy surfaces labeled by a time 
function t, and let t"" be a time evolution vector field, satisfying f^Vat = 1, which may be 
decomposed into a lapse and shift via 

e = Nn" + N" , (2.8) 

where n"" denotes the unit normal to T,t, and where A^" is tangential to Sj. As noted above, for 
any $ S W which is C^, the quantity {n"'Va^)\T:t is uniquely determined by ^{sf Assuming 
that (n"Va^)|st hes in L^(Si) for at least a dense (in L^(St)) subspace of <1> G D, we obtain, 
for each t an operator Ct : L^(Sj) -^ L'^(T,t) such that 

(n'^V„<i>)|E, = -ia(cl>|sJ. (2.9) 

In the following we shall omit writing "Ist", i.e., it should be understood that, for example, 
Ct$ means Ct($|x;j). It follows immediately from Eqs. ( p.4| ) and ( |2.5| ) that for <I>i,<&2 in the 
subspace of T> on which Cj is defined, we have 

(^l,^2)/rG = i(^l,n"Va^2)L2{St)-«(^''Va$l,$2)L2{S0 
= ($1, Ct$2)L2{E0 + (Ct^l, $2)l2(E0 

= i<^i,{Ct + Cl)^2)m^,). (2.10) 

By comparison with Eq. ( p.6D , we obtain 

Bt = {Ct + Cl)'/\ (2.11) 

By brute force substitution, we obtain as the evolution equation for ^t = Bt^t 



dt dt 



{Bt<^ 



'^^-^.^B.'^^ 



dt ' ' dt 

B^^^t - iBtNCtBi'^t + BtN'^VaiBi^^t) 



dBt 



dt 
-iHeff^t, (2.12) 



where the effective Hamiltonian for evolution in L^(St) is given by 



dB 

Heff = i-QfBi^ + BtNCtBi' + iBtN^VaBi^ . (2.13) 



Since Bt is given in terms of Ct by Eq. (2.11), the effective Hamiltonian is determined directly 
by Ct- Note that since the natural volume element on Sj will, in general, be "time dependent" 
(i.e., not Lie-derived by f"), the natural L^ inner product on Sj will similarly be time dependent, 
and H(,ff will not be self-adjoint. However, we could define a self-adjoint evolution by working 
with a fixed "coordinate" volume element on Sj to define a fixed inner product for L^(Ei) or, 
equivalently, working with densitized versions of ^f We have chosen not to do so here since 
this would further complicate the formulas, and, in all of the applications in this paper, the 
natural volume element on S^ will be Lie derived by our time evolution vector field i". 

The above general relations simplify considerably for the quantum cosmological models we 
shall consider in this paper. First, we shall choose as our Cauchy surfaces which foliate mini- 
superspace the hypersurfaces orthogonal to a timelike Killing field t". We choose t"" as the time 
evolution vector field, so we have A^" = and the Cauchy surfaces are labeled by the time 
function t such that f^Wat = 1. Furthermore, in our models, we have A^ = 1. Most importantly, 
in our models the f-dependence is separable, so we may work with (un-normalizable) basis 
functions for TC of the form 

,^^^(t,q^) = ^S^M), (2.14) 



where fuj{t) is such that fuj{t) -^ g-i^t ^s t ^ — oo. Here g* are the configuration variables apart 
from t, and the label a represents the quantum numbers other than uj. The normalization of 
ScuaiQ^) is chosen to be (5^/o-', S^a)L'^ = (^(f^' — io)5{a' , a), where 6{a' , a) is an appropriate delta- 
function for the label a, so that (^oj'ct') ^u)a)KG = ^{yj' — uj)6{a' , a). Note that the conservation 
of the Klein-Gordon inner product together with the asymptotic form of fuj{t) as t — > — oo 
implies that 

Ut)fUt) - Ut)fL{t) = -2^w , (2.15) 

where /' = df/dt. 

Differentiating Eq. ( p. 14 ) with respect to t, we obtain 

d^u,. fUt) 



dt y/2^ 

fUt) 



^ua ■ (2.16) 



Ut) 

Thus, by inspection, the operator Ct defined by Eq. ( |2.9D is diagonal in this basis, and we have 



Using Eq. ( |2.15|) , we obtain 



Ct^u^a = ij^'^ua . (2.17) 



Bt<^ua = (Ci + C/)i/2cD, 

y/2uj 



lutf- ''^''' 



Substitution of the above expressions for Ct and Bt together with A^ = 1 and N'^ = into Eq. 
( 2.13| ) then yields 



H,ft^ 



en ^ wa 



<^ ( 1 \,, „M , J'M 



•s(^)i-f"("i+' 






ijj 



\Ut)\ 



Ut). 

■^u^a , (2.19) 



where Eq. (2.15) was used again in the last line. 



3 Comparison with Deparametrization 

An alternative scheme to the one presented in the previous section for constructing a quantum 
theory corresponding to a classical theory with super-Hamiltonian of the form ( |2.1| ) is to "de- 
parametrize" the theory prior to quantization in the manner proposed by Arnowitt, Deser, and 
Misner [g]. In this procedure, we first choose a time function, t, whose level surfaces are Cauchy 
surfaces, Sj, and choose a time evolution vector field, t", satisfying t^'Vat = 1. We then solve 
Eq. ( |2.lD for the momentum, ttj, canonically conjugate to t, thereby expressing nt as a function 
of t and the remaining dynamical variables (g*, ttj). The Hilbert space of states at time t is 
then taken to be -L^(S(), and the quantum operators corresponding to g* and ttj at time t are 
defined by the standard prescription for L^(S(). Dynamical evolution then is defined by 

— = -iHADM^, (3.1) 

where the ADM Hamiltonian is given by 

Hadm = -vrt(t, q\ ttj) . (3.2) 

Classically, this Hamiltonian generates the dynamics of {g'*,7rj}, with respect to time t, as can 
be verified by using the original Hamilton's equations with the super-Hamiltonian H. 

Since the proposal of the previous section also can be viewed as a theory defined on the 
Hilbert spaces L^(S() with the standard definitions of position and momentum operators, it 
differs from the quantum theory obtained by the above deparametrization method only via 
the use of the Hamiltonian (2.13|) rather than (p.2|). In this section, we shall compare the two 



approaches for the "toy model" of a relativistic particle in four dimensions with a potential of 
compact spacetime support, whose super-Hamiltonian is given by 

H = -pI + ^'^ + \V{xq,^) = Q. (3.3) 

We study this model to first order in a perturbation series in A. We shall show explicitly 
in this model that - as expected (see, e.g., [^) - the dynamical evolution defined by the de- 
parametrization procedure depends upon the choice of time variable, t, used to perform the 
deparametrization. The quantum theory defined by the proposal of the previous section is 
manifestly independent of a choice of slicing, but we shall show explicitly the presence of a 
"memory effect" , which illustrates the nonlocal-in-time character of dynamical evolution. 

First we solve the Klein-Gordon equation ( |2.3| ) for this model to lowest order in A. We have 
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^ + ^h> = Anxo,x)<.. (3.4) 



Consider the "unperturbed" wavefunction 



^^°^=/(2^^(P)^'^'"^^'"'' ^^-'^ 



where p = |p|, and write the sohition to ( p. 41 ) to first order in A as 

$ = $(0) + $(1) (3.6) 

with the boundary condition ^^^' -^ for xq -^ — oo. Note that the Klein-Gordon norm of ^^^' 
is 

The first order in A correction, $(^-', to ^ due to the presence of the potential is given by 
$«(xo,x) = Ay'd4dVG^G(a;o,x;x'o,x')yK,x')$(°)(x'o,x'), (3.8) 

where Gkg denotes the retarded Green function for the Klein-Gordon operator, 



Gkg{xo, x; Xq, x') = -ie{xo - x'q) J ^^^^ ^ 



ip' {xo-x[^)+ip' ■{:x.~x.') _ ip'(xo-a:Q)+jp'-{x-x') 



g 'f \-u -„y . -f \- ^ ' — e 



(3.9) 

Let Xq denote the maximum value of xq such that V{xq,x.) ^ for some x. Then we have 
for xo > Xg , 

cl>«(xo,x) = Ay'(ix^dVGxG(xo,x;x'o,x')F(x^,x')$(°)(x'o,x') 

-^' 1 1 (2;)32p(2;)32p4 ^(^^-^'p-p)^(p)-"""°"^^-^ 

-V{-p' -p,p'- p)/(p)e*^'-«+*P'-l , (3.10) 



where 



V{po, p)= f dxod^xV{xo, x)e^P°^°-^P-^ . (3.11) 



Consider now the quantum theory obtained by deparametrizing with respect to the time 
variable t = xo- Then we have 

7rf =p'^ + XV {t,:xi). (3.12) 

The first point to note is that if V becomes negative, a real solution, ttj, of Eq. (3.12) will not 
exist for all values of t, x and p, and the quantization prescription will break down. In order to 
avoid this difficulty, we restrict consideration to potentials which are everywhere non-negative. 
(Note that, in contrast, the proposal given in the previous section does not require any such 
restriction on V.) If V is non- negative, we have 



Hadm = -vTi = Vp' + Al^(t, x) . (3.13) 

One of the nice features of this model is that the definition of Hadm as a quantum operator is 
unambiguous: no "factor ordering" ambiguities occur for defining the operator p"^ + XV, and, by 
continuity, one must choose the positive square root of the operator in order to get agreement 



with the standard free particle theory at A = 0. Thus, the dynamical evolution equation in this 
approach is 



■ ^ T 

t^ — "^ADM 
OXq 



52 



1/2 



^ ADM ■ 



Now consider the solution to Eq. (3.14) to first order in A of the form 



^ADM = "^foM + *aJ)M 



with \1' 



(1) 
ADM 



for xq -^ — cxD, where 



^ 



(0) 
ADM 



(fl 



(27r)3V2p 



/(P)e 



-jpxo+ip-x 



(3.14) 



(3.15) 



(3.16) 



Note that '^j^jj]\^ corresponds precisely to ^^^' , Eq. ( |3.5| ); the difference in normalization arises 
from the fact that "^y^^j^j^f is normalized via the L^ inner product, whereas ^^^' was normalized 
via the Klein-Gordon inner product. 

The square root operator in ( 3.14| ) acts on e*P'-^ to order A as 






1/2 



gipx ^ ipx 



pe^P-^ + A 



(27r)3 p + p' 



,«P x 



where 



y,(a;o,p)= /d3xy(xo,x)e-^P-'' 



(3.17) 



(3.18) 



One can readily verify that the square of this operator is indeed —d'^/d'K^ + Ay(a;o, x) + 0{}?). 
From Eqs. ( 3.14 ) and ( 3.17 ), we find that the first order correction, ^ j^£){^t to ^ adm due to 



the presence of V satisfies 



d 



' dxo 



d-x? 



1/2N 



^^M = A / ^'P ^-^"°'P\-PV (p)e-P'- . (3.19) 



Using the retarded Green function for the Schrodinger operator appearing on the left side of 
this equation 



G^dm(2;o,x;xo,x') = -ie{xQ - x'q) 



d' 



(27r)3 



P_ -jp(a;o-x^)+jp-(x-x') 



(3.20) 



we find for xq > Xq 



* 



(1) 

ADM 



-iX 



(fp d^p' V{p' — p, p' — p) 
(27r)3y2^(27r)3 p + p' 



fip)e 



-ip a;o+2p -x 



(3.21) 



In the region xq > Xq , we can associate to '^ adm a positive frequency solution '^adm to 
the free Klein-Gordon equation by the correspondence 



^ADM = (-4V2)-V4^^^^ . 



(3.22) 



Writing '^adm = ^adm + ^'^adm + ^(A^), we have 



^^DM = '^^°^ (3-23) 

(27r)3^/2^(27r)3y2^ p + p' 



^%M = -^A / ^^:j^,^:j^''' :^'^, -^v (p)e-^'-o+-p'-, (3.24) 



where <I>(°) was given by Eq. ( |3l5|) . Note that ^ajdm differs from <^(^) in Eq. (|3.10|) . Indeed 



unhke '^ ajj^^ ^ ^^ ^°^ even a purely positive frequency solution. However, we can consider 
a WKB expansion of ^a£>m ^'^'^ ^^^ ■ For V smooth, the negative frequency part of <^'^^ 
is "nonperturbative" in such an expansion, i.e., it does not appear to any finite order in fi. 
Furthermore, we have verified that to first order in h, ^adm ^'^'^ ^ agree. This result is 



consistent with Barvinsky's [10| much more general arguments in the context of quantum field 
theory that "reduced phase space quantization" agrees with "Dirac quantization" at one-loop 
order. However, in our model, we find that <&^£)^ and ^'^-^ differ at order Ti . 

Now we consider a second deparametrization scheme obtained by using a boosted Lorentz 
frame, i.e., we choose t = x^ = xocosh/? + xisinh/3. The initial ADM wavefunction which 



,(0) 



corresponds to ^ aj:,]^ given by ( 3.16 ) is 



"^ADM-J ^^^^s^nP~(s)e , (3.25) 

where x^ = (xi cosh /3 + xosinh/3, 0:2, ^3) and p_/3 = {picoshp — posmh.p,p2,P3) ■ If we start 
from this initial wavefunction, we obtain for the lowest order correction, "^ ADJ^^, an expression of 



the form ( 3.21| ) with xq replaced by Xq = xo cosh /? + xi sinh /3, xi by xf = xi cosh P + xq sinh /3, 



/(p) by fip-fs) and V{p' -p, p' - p) by F(p'_^ -p-p, p'_/3 - P-fs)- The corresponding solution, 
^ADM^ of the free Klein-Gordon equation at late times is given by 



(UP r d^p (fp' '^Jppp'pVip' -p,p' -p) ■ , ^- , 

^^^^- '-^y (2vr)32p(2vr)32p' pp+p'^ ^^P^' ' ^^'^^^ 

where pp = {pi cos\il3+pQsmh.f3,p2,Pz) andp^ = \pf3\. This disagrees with the solution '^a£im, 



Eq. ( |3.24 ), obtained with the original choice of time slicing. Hence, if we physically identify ADM 



wavefunctions on different slices in the region to the future of the potential if they correspond 
to the same Klein-Gordon wavefunction, we find that the final state obtained here is physically 
distinct from the one obtained with the original choice of time slicing. Thus, this example — 
which is not clouded by factor ordering ambiguities - explicitly shows the dependence of the 
deparametrization method upon the choice of time slicing. It is also possible to show slice 
dependence of the deparametrization method with fixed initial and final time slices in lowest 
order in perturbation in the deformation of slices with a natural factor ordering [^] . 

Next, we study the model ( p.3[ ) using the quantization procedure outlined in Sec. |2|. We 
shall study the dynamics predicted by this procedure by expressing the theory in L^ form as 
described at the end of Sec. |2|, using the time slicing t = xq, so that a direct comparison can be 
made with the deparametrization procedure. As we shall see, a "memory effect" is present: The 
dynamics defined by Hf>^ff does not return to standard, free particle dynamics after the potential 
has been turned off, i.e., the system has the "memory" that it had a nonzero potential in the 
past. Although this effect will not be relevant for the cosmological models that we will study in 
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this paper, it could be important, for example, in models in which the universe tunnels through 
a potential barrier in an early stage. 

In order to calculate H^jf, we must calculate the operator Ct defined by Eq. ( ^^ ) above. 



S^) 



Since by Eq. ( 3.10| ) the Klein-Gordon solutions in the region xq > Xq"'' which are asymptotically 
ositive frequency in 



positive frequency in the past take the form (to first order in A) 



{2T:f2p 



y{p - P, P - p)e-^P'*+*P'-^ - V{-p' -p,p - p)e^P'*+^P'-'' 



we find for xq> Xq 



Cte'P'' = pe'P-'' - i\ 



^y(-p'-p,p'-p)e^(^+P')*+^P'-. 



(3.27) 
(3.28) 



It follows that to first order in A the operator Bt = [Ct + CJ]^'"^ is given by 



Bte'P'' 



2pe 



ipx 



+iX 



d^p' V{p' +p,p' - p)e-*(P+P')* - V{-p' -p,p' - p)e'(P+P')* ^p,., 
(2^ '' 



2p + V2F 



(3.29) 



Since A'^ = 1 and A^" = 0, for xq > Xg we obtain from Eq. ( |2.13| ) the result 

TT „ipX „„JpX 

^effe^ — pe^ 

+iX 



d^p' 



(2vr) 



Vp 



y/P+VP 



=V{p' +p,p' - p)e-*(p'+P)*+*P'-^ 

i{p'+p)t+ip'-x 



Vp + 



V{-p'-p,p'-p) 
yV 



(3.30) 



It is readily apparent that Hf.f^ ^ Hadm- Furthermore, although Hejj agrees with the standard, 
free particle Hamiltonian Hq = V— V^ prior to time when the potential is "turned on" , we have 
Heff 7^ Hf) after the potential is "turned ofF'. This means that, in principle, by observing the 
dynamics of "free" particles (after the potential has been turned off), one could deduce that 
a potential had previously been present. This "memory" phenomenon can occur because the 
positive frequency condition in the asymptotic past enters the definition of the operator Ct - 
and, thereby, H^jf ~ thus making these operators nonlocal in time. 

Some insight into the nature of this memory effect can be obtained by considering the simpler 
case where the potential is purely a function of time, t (and, thus, is not of compact spacetime 
support). Then the positive frequency solution which in the past (before the potential is turned 
on) behaves as g-^pt+^p-^ evolves in the future (after the potential is turned off) to a solution 
of the form 



-jpt+*p-x 



a 



^p)g-ipt+ip-x ^ ^(p^g+ipi+ip.x ^ 



(3.31) 



JA/) 



where |a(p)|^ — |/9(p)p = 1. Then, it can be readily seen that for xq > Xq"' the effective 
Hamiltonian is 



H^ffe^^-- 



p 



|a(p)e-^P*+/3(p)e^P*|2 



=«PX 



(3.32) 
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i.e., after the potential has been turned off, the Lp' wavefunctions evolve as 

|a(p)e *P* + p(p)e*P'| 

In this case, the phase acquired in the period T = 2t:/p is still — 27r since |a(p)| > |/3(p)|- This 
implies that for a wavepacket sharply peaked around p, the motion of the expectation value of 
the position operator, x, will simply oscillate about the standard free-particle value. However, 
more complicated dynamical behavior would occur in the model considered above where V has 
compact spacetime support. 

4 Homogeneous cosmologies 

In this section, we will provide more details concerning the application of the proposal outlined 
in Sec. |2| to general homogeneous cosmological models with a scalar field [^]. The possible 
homogeneous cosmologies are comprised by the Bianchi models and the Kantowski-Sachs model 
(see, e.g., |1^). We first consider the Bianchi models, following the discussion and notation of 
Ref. |1J]. In the Bianchi models, the spacetime manifold, M, is taken to be IR x G, where G is a 
three dimensional Lie group. (We parametrize IR by the variable, t.) For any k € G, we define 
a diffeomorphism ipk - called left translation by fe - by ipk{g) = kg for all g £ G. A tensor field 
T"""b... is said to be left invariant if ip%T°""b... = T"'"'b..., where Vfc is the map on tensor fields 
induced by ipk- The spacetime metric on M is taken to be of the form 

gab = -yatVbt + habit), (4.1) 

where the spatial metric, hab, on G is left invariant, and thus can be expanded as 

3 

hab=J2^^j(t){a'Ua^)b, (4.2) 

where (cr*)(j are left invariant one-forms that are independent of t. 

The structure constant tensor field C"^^^ on G is defined by [v, wY = C'^^^^v°"ufi for any two 



left invariant vector fields v"" and w"" . It is known that G'^^^^ can be expressed as |15| 

C\, = M-''eaab + 5laAb], (4.3) 

where M is a left invariant symmetric tensor on G, Aa is a left invariant one-form and eabc is a 
left invariant three-form on S. The Jacobi identity implies that M°'^Ab = 0. The Lie algebra is 
said to be of class A.\i Aa = 0, and of class B otherwise. No consistent Hamiltonian formulation 
is known for the Bianchi models of class B |16|. In particular, the Hamiltonian constraint does 



not generate the dynamics of the system. We specialize to the models of class A for this reason. 
The Lie algebras of class A are uniquely determined up to isomorphisms by the rank and 
the signature of the tensor M , where the overall sign of M is irrelevant. Thus, one has the 
following six inequivalent cases: 

(0, 0, 0), (+, 0, 0), (+, -, 0), (+, +, 0), (+, +, -), (+, +, +). 

The corresponding Lie algebras are called the Bianchi types I, II, VIq, VIIq, VIII and IX, 
respectively. The unique connected, simply connected Lie group G for Bianchi type IX is 
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SU{2), which has topology S^. The corresponding simply connected Lie groups for the other 
Bianchi types are nonconipact, and only the (trivial) Bianchi I algebra is also the Lie algebra 
of a compact Lie group (namely C/(l) x C/(l) x U{1)) [ p^ . 

In the compact case, we choose the configuration variable, hab, for the gravitational degrees 
of freedom to be the spatial metric at the identity element, e. The conjugate momentum 
variable vr'^*, is then the usual momentum density at e multiplied by the "volume of space", 
V, as determined by the 3-form Cabc introduced in eq.(O) (so that vr'** may be viewed as the 
integral of the momentum density over space) . Similarly, we choose the configuration variable 
for the scalar field to be the value of (/) at the identity element and its conjugate momentum, vr^, 
is taken to be the usual momentum density multiplied by V. The momentum constraint reads 

Pa = Tr'',M^'edab = 0. (4.4) 

The Hamiltonian constraint is |15] 



(4.5) 
The factors for the scalar contribution have been chosen for later convenience. Equations ( |4.4D 
and ( |4.5| ) also hold in the noncompact case (with V an arbitrary constant) after suitable (infinite) 
rescalings have been made on the dynamical variables and H. 

In Appendix A, we shall consider the quantum theory obtained by treating the constraints 



(4.4) and (|4.5| ) on an equal footing. In particular, we shall derive the explicit form of the 



Wheeler-DeWitt equation for the Bianchi IX model which results from imposing (4.4) as an 
operator constraint on the wavefunction. However, in the body of this paper we shall consider 
the models obtained by first reducing the problem at the classical level - thereby eliminating 



the momentum constraints classically - and then imposing ( [4. 51) as an operator constraint. We 



achieve this reduction by noting that the classical momentum constraint (4^) is equivalent to 
the statement that at each time, t, the quantities ir'^b and M'^f, commute when viewed as linear 
maps acting on the tangent space, Ve, at the identity element e £ G. (Here, all indices are 
raised and lowered with hab and /i'*''.) This means that we can simultaneously diagonalize these 
linear maps in an orthonormal basis of hab- The classical evolution equations then imply that 
7r°*, M"^, and hab will remain diagonal in this basis for all time. Thus, classically, there is no 
loss of generality in restricting to the "diagonal case" , in which case the momentum constraints 
become trivial. 

It is convenient to choose the diagonal basis ((T*)a so that Mi = M"'^{a^)a{(J^)b takes the 
form 

(Mi,M2,M3) = (0,0,0), (+1,0,0), (+1,-1,0), 

(+1,+1,0),(+1,+1,-1),(+1,+1,+1), 

respectively, in the six cases corresponding to the different possible choices of signature. We 
then write the basis components of hab and ir""^ as hij = diag(/ii, /12, /13), vr*-' = diag(7r^,7r^, tt'^) 
and hi = h~^'^hi. We define a, /?+, and /3_ by 

hi = e2a+2{/3++v/3/3_) ^ (4^g) 

h2 = e2"+2(/5+-v^/3-) , (4.7) 

/j3 = e2°-4^+. (4.8) 
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(Note that this definition of a is consistent with Eq. (2.2), i.e., the maps on niinisuperspace 
defined by Eq. ( |2.2D correspond to a translation of the parameter a defined by the above 
equations.) In terms of these variables, the Hamiltonian constraint becomes 



B. 



V 



1/2 



e-^"(E: + p) = 0, 



where 



K 
P 



-vTq, + tt^ + 7r_ + vr j, 



(4.9) 



(4.10) 



Aa 



Mfhj + M^hl + M|/i; 



-M2MJi^^ - M^mJi^'^ - MiM2h^^ 



+e*'"^.(</') 



(4.11) 



where tTq, it± and vr^ are the conjugate momenta of a, f3± and (p, respectively. We also have 
eliminated the factor of 24V^ in P by shifting a by (1/2) ln(\/24V). The Wheeler-DeWitt 
equation corresponding to (|4.9| ) then is simply a Klein-Gordon equation in flat spacetime with 
potential P, i.e., 



di 



dl 



dl 



di + P 



$ = 0. 



(4.12) 



In the Kantowski-Sachs model ||l^, space has the topology of S^ x S"^ acted upon by the 
isometry group U{1) x 50(3). After performing a similar, classical reduction of the momentum 
constraints, the spacetime metric of this model can be written in the form 



ds^ 



.df + e2["W+2/5W]rf^2 ^ ^2[a(t)-m] (^^2 ^ ^-^2 



(4.13) 



With scalar field matter present, the rescaled super-Hamiltonian takes the form (after some 
redefinition of variables) 



H 



-7r„ + TT^a + vr^ 



Ja+2fS^^6ay 



0, 



so the Wheeler-DeWitt equation is simply 

5/3 -f^,/) 



dl-dl-dl + P 



$ = 0, 



4a- 



-^f + e'"Vs^ 



(4.14) 



(4.15) 
(4.16) 



It is worth elucidating the origin of the differences occurring between the quantum theory 
obtained by the above classical reduction of the momentum constraints as compared with the 
method described in the Appendix A. For definiteness, we focus attention upon the cases of 
the vacuum Bianchi I and Bianchi IX models. In the Bianchi I case, we have M = 0, so the 
momentum constraint is absent entirely. The above "classical reduction" to the diagonal case 
is closely analogous to formulating a quantum theory of a free, nonrelativistic particle moving 
in IR by first noting that, classically, the particle moves in a straight line, so that, classically, 
we may choose our coordinates in IR so that the particle moves only in the x-direction. We 
may then formulate a quantum theory equivalent to that of a free particle moving in IR . 
However, this reduced theory - analogous to the treatment of the Bianchi I model given above 
- is not equivalent to standard free particle theory in IR , which is the analog of the treatment 
of the Bianchi I model given in Appendix A. On the other hand, for the Bianchi IX case, the 
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momentum constraints are nontrivial, and the number of degrees of freedom are the same for 
the above classicahy reduced theory as for the theory obtained via a quantum reduction of 
the momentum constraints as given in Appendix A. A good analog of the Bianchi IX model is 
a free, nonrelativistic particle in IR with the imposition of the additional constraint that its 
vector angular momentum vanishes. The analog of above classical reduction procedure would 
be to note that, classically, the particle moves on a straight line through the origin, and then 
to formulate a quantum theory equivalent to that of a free particle moving in IR . The analog 
of the procedure given in Appendix A would be to solve the constraint by restricting attention 
to wavefunctions in IR which are spherically symmetric. This yields a theory which is not 
unitarily equivalent to standard free particle theory in IR |20]. 



Clearly, the quantum treatment of the momentum constraints given in Appendix A is "more 
correct" than the classical reduction of them given above. Nevertheless, in the body of this 
paper we choose to work with the models obtained by the above classical reduction, since the 
Wheeler-DeWitt equation is considerably simpler in this case - being a wave equation in a flat 
rather than curved spacetime. In addition, most of the previous literature has considered the 
Bianchi models obtained by the classical reduction, so our consideration of them here should 
facilitate comparison with previous approaches. Since the models will be used only to elucidate 
the qualitative features of the quantum theory proposed in Sec. § - rather than to attempt to 
make realistic predictions about the actual universe - we see little disadvantage to working with 
the simpler, classically reduced models. 

We turn our attention, now, to an investigation of whether the following two conditions - 
which are necessary and sufficient for the implementation of the proposal of Sec. y - are likely 
to be satisfied in the above Bianchi and Kantowski-Sachs models: (1) For the construction of 
the Hilbert space, Ti., of states, it is necessary for the potential term in the Wheeler-DeWitt 
equation to vanish sufficiently rapidly as a — > — oo that the notion of "asymptotically positive 
frequency solutions" is well defined. (2) For the construction of observables on a "time slice" 
Sj, it is necessary that the space, V, of C^ solutions in Ti whose restriction to S^ lies in L'^(T,t) 
be dense in both TC and L^(Sf); it also is necessary that any sequence in V which converges in 
both TC and L^(S() have a nonzero limit in TC if and only if it has a nonzero limit in L^(S() 0. 

With regard to condition (1), it should be noted that in all of our models (as well as in full 
quantum gravity), the potential terms vanish exponentially rapidly as a ^ — cxd. Nevertheless, 
since the potential terms also may blow up exponentially rapidly at spatial infinity, it does not 
automatically follow that these terms can be neglected at "early times", in the manner needed 
for the construction of TC. 

The following criterion is very useful to consider for the analysis of condition (1): Consider 
the vector space, S, of real solutions to the Wheeler-DeWitt equation with initial data of 
compact support on Cauchy surfaces. Define an "energy inner product" on S associated with 
each Cauchy surface a = const, by 



J a=const. 



dA^ldB<^2 + \GAB{dc<^id^<^2 + P'^1'^2] 



\B 



[d/day , (4.17) 



where Gab is the flat metric appearing in ( f4.12 ) and ( [4.15| ), so that the "energy norm" at "time" 



a of a solution, $, is just the integral over the surface at that value of a of the stress-energy 
tensor of <1> contracted once with the unit normal and once with d/da. Then the following 
condition is sufficient to enable the desired Hilbert space, TC, to be constructed: 
Condition (1' ): For all <^i,<l>2 G S the limit 



E($i,$2)= lim E^{^i,^2) (4.18) 
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exists and satisfies 

E($i, $i)S($2, $2) > K\n{<^i, $2)1^ (4.19) 

for some constant K > 0, wliere Q was defined by Eq. (p.4|). 

That condition (1') implies tliat tlie desired 7i can be constructed follows immediately from 
the work of Chmielowski [pi]]. We define the inner product, fi, on S satisfying 

M<I>i,<5i)- *.^0 4^(^^^^^) (4.20) 



to be the inner product associated to E by the construction of Proposition 1 of pTI]. A Hilbert 



space, Tl, of solutions then can be constructed from // as discussed in detail in, e.g., |]2^ and |23]. 
This Hilbert space has the interpretation of being comprised of the "asymptotically positive fre- 
quency solutions" for the same reason that the standard one-particle Hilbert space for stationary 
spacetimes [^ has the interpretation of being comprised of positive frequency solutions. 

In the vacuum case, we believe that condition (1') holds for the Bianchi models I, H, VIq, 
and VHq as well as for the Kantowski-Sachs model. However, it appears that condition (!') 
will fail for the vacuum Bianchi VHI and IX models. Nevertheless, we believe that condition 
(1') will hold in all of the homogeneous cosmological models (including VIII and IX) in which a 
scalar field, (j), is present, provided only that the potential, V{(j)), does not grow exponentially 
rapidly in (p. 

Our arguments for these beliefs are based mainly upon the classical "particle" dynamics 
associated with the super-Hamiltonian ( [4. 91) , since the behavior of suitable "wavepackets" satis- 



fying ( 4.12| ) should be similar to this classical particle dynamics in the limit a ^ — 00. Consider, 



first, the vacuum case. The spatial region in the (/?+—/?_ )-plane where the terms 



„4o/ n,f_!\,f_I-l (i/T^d/T. ?^-l ]\J. ]\/r^h-'^'\ 



e*"(-M2M3/i7' - MsMih^' - MiM2/i3 j 



in P of ( [4.11 ) are large "moves away" from the origin at the speed of light as a ^ —00. Hence, 



for large enough —a, these terms should be negligible in both the particle and wave dynamics, 
and the relevant contribution from P should be 



1 



P « -{Mfhl + Mfhl + Mfhl 



^ :^4a+4/3+ +4^3/3. , ^ 4q+4/3+-4V3/3_ , ^ 4a-8/3+ (^ 2I) 

2 2 2 ■ \ ■ J 

These terms are well approximated as a ^ — oo by potential walls located, respectively, at 
/?+ + \/3/3_ = — a, /3+ — \/3/3_ = —a and —2/?+ = —a [24|. These walls surround a region in 
the shape of a triangle and recede from the origin at half the speed of light. In Bianchi types I, 
II, VIq, and VIIq, the rank of M"'^ is less than three, i.e., one or more of the Mj's are zero, so 
at least one of these walls will be missing. In that case, if we consider evolution backwards in 
a, a generic particle trajectory should bounce from the remaining walls at most a finite number 
of times, after which time it can be treated as a free particle. Similarly, wavepackets should, 
after a finite time, escape to a region where they satisfy the free Klein-Gordon equation to an 
excellent approximation. Condition (1') should then be satisfied in the vacuum Bianchi I, II, 
VIq, and VIIq models. (Indeed, these arguments can be made rigorous for the vacuum Bianchi 
I and II models, since they can be solved explicitly.) Similar arguments apply (rigorously) to 
the vacuum Kantowski-Sachs model, where the region where the potential, — e"^"^^^, is large 
recedes from the origin at twice the speed of light. 
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In Bianchi type VIII and IX models we have Mf = M2 = M3 = 1, so all the "potential 
walls" are present. Since these walls recede only at half the speed of light ~ whereas in the 
classical dynamics, the "particle" in the minisuperspace always moves at the speed of light 



between "bounces" - an infinite number of "bounces" generically occurs as a — > — cxd |25]. 
Since energy is lost on each of the "bounces" the energy should asymptotically approach zero. 
Similarly, in the wave dynamics, it appears that Eq. — > as a ^ —00, and condition (!') should 
fail to hold. 

However, the situation improves considerably if scalar field matter is included, provided that 
V{(j)) does not grow exponentially rapidly in \(j)\. In that case, the scalar potential term, e^^V^cp), 
should become negligible for large —a, and the momentum vr^ should become approximately 
conserved. The classical dynamics as a — > —00 then corresponds to that of a massive particle, 
with mass 7r|. Although the potential walls still recede at only half the speed of light, on 
account of the mass term the classical particle now will slow down each time it bounces off 
a wall. Eventually, its velocity will drop below half the speed of light, at which point the 
particle will no longer see the potential walls. Thus, if vr^ 7^ 0, only a finite number of bounces 
will occur, and the classical dynamics will not be chaotic as a ^ —00, i.e., the presence of a 
scalar field is sufficient to eliminate the classical chaotic behavior of the Bianchi VIII and IX 



models [26|. Furthermore, the energy of classical solutions will be approximately conserved as 
a — > —00. Similar behavior should occur for wavepackets, so it appears that condition (!') 
should hold. Thus, it seems likely that all of the homogeneous cosmological models with a 
scalar field will satisfy the first condition needed for the implementation of the proposal of Sec. 
g. Note, however, that it is known that the classical chaotic behavior is not eliminated if one 
includes a homogeneous electromagnetic field rather than a scalar field in Bianchi type IX model 



1 27 1, so it appears that a scalar field is a necessary ingredient for the construction of the Hilbert 
space in the Bianchi VIII and IX models. 

As for condition (2) - which requires that there be a dense subset of L^-normalizable func- 
tions in the Hilbert space - we see no reason to doubt its validity in the Bianchi models for 
the surfaces of constant a. However, in the Kantowski-Sachs model the solutions to Eq. ( [4.15| ) 
in the vacuum case which can be obtained by separation of the variables r = (2a + /3)/v3, 
^ = (2/5 + a)/\/3 turn out to grow faster than exponentially for large /5. We expect condition 
(2) to hold on the surfaces of constant r, but the rapid growth of the solutions in f5 suggests that 
condition (2) will fail for the surfaces of constant a, although we have not been able to prove 
this. Similarly, in the Bianchi IX model (as well as models obtained by suppressing some of its 
degrees of freedom, such as those studied in the next two sections), it appears that condition 
(2) will fail for some choices of Cauchy surface other than the surfaces of constant a. 

5 Robertson- Walker universe with a scalar field 

In this section we apply the quantization prescription of Sec. § to the closed Robertson- Walker 
cosmology with a homogeneous, free, massless scalar field. We shall choose the hypersurfaces 
of constant a as our time slices in minisuperspace. 

The classical (rescaled) super-Hamiltonian can be obtained from that of the Bianchi IX 
model, Eq. (^), by setting tt± = P± = 0, i.e., 

^W = -vr^ + vr^-e^" = 0. (5.1) 

The trajectories in minisuperspace corresponding to the classical solutions of the equations of 
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motion for this super-Hamiltonian are given by 

a = i {In \p\ - In cosh [2(,/. - </.o)]} , (5.2) 

where (po and p (= tt^) are arbitrary real constants. These trajectories describe classical uni- 
verses which start at a "big bang" singularity at a — > — oo (with \(p\ -^ oo), expand to the 
maximum size om = (1/2) In \p\ at (j) = (po^ and then recollapse in a symmetrical manner. 

We have two main motivations for studying this model: (i) As noted above, all of the 
classical solutions recollapse. On the other hand, the proposal of Sec. y treats a as a "time 
variable", which can be prescribed arbitrarily to "set the conditions" for the other dynamical 
variables. Thus, in particular, we could choose a state in 7i which, at a ^ — cxd, corresponds 
closely to a classical trajectory (^^), and ask about the behavior of for a >> om- How does cj) 
behave in this classically forbidden region of minisuperspace? (ii) The classical trajectory (|5.2D 
is everywhere spacelike in minisuperspace. On the other hand, the propagation defined by the 
Wheeler-DeWitt equation is entirely causal with respect to the light cones defined by the metric 
on minisuperspace. Thus, one might anticipate some difficulties in approximating a classical 
trajectory by a state in TC even during the "expanding phase" of the classical solution. Clearly, 
the states in TC do not behave classically for a > om- How close can one get to om before the 
universe begins to behave nonclassically? 

As we now shall see, for wavepackets in TC corresponding to classical solutions which expand 
for much longer than the Planck time, the answer to question (i) is that for a > um , 4' "freezes" 
at the value, (po, corresponding to the classical value of (j) at maximum expansion, a = um- The 
answer to question (ii) is that despite the spacelike character of the classical trajectory, suitably 
chosen wavepackets in TC accurately describe the expanding phase of the classical solution until 
a K, aM- 

The Wheeler-DeWitt equation for this model can be obtained by setting vr^ = —ida and 



TT^ = —id(f) in Eq. (5J.), i.e., we have 



Hrw^ 



da^ 



g4a 



$(«,(/>) = 0. (5.3) 



Note that the deparametrization approach of Sec. ^ cannot be employed here because the 
potential is negative, but the method of Sec. |2| can be used without difficulty. 

Since the momentum vr^ is conserved, the solutions to the equation Hrw^ = can be 
written as $p(a,(/)) oc /|p|(a)e*^'^, where the function f\p\{a) satisfies 




„4« 



e 



/|p|(a)=0. (5.4) 



The two linearly independent solutions of this equation are /_|_j|p|/2(e2'^/2), where 

M-)-E „r(. + . + i) (2) ■ f^'^) 

The asymptotic positive frequency condition requires that /|p|(a) -^ e"*!^!" for a — > —00. Thus, 
the solutions relevant for constructing wavepackets which lie in TC are 

/|,|(a) = 2-IPlr(l - z|p|/2)/_,|,|/2(e2V2) . (5.6) 
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If we impose the (5-function normalization {^p,^p')KG = ^{p ~p')-, our "basis" for solutions in 
His 



The function Iu{z) behaves for large \z\ as |2^ 



h{z) ~ -= . (5.8) 



e 

llTZ 

Hence, in the classically forbidden region e" >> \p\, the function f\p\{a) grows like exp(e^°/2 - 
a). This behavior of /|p|(a) can be understood as follows. Eq. (|5.4D has the form of a time' 
independent Schrodinger equation, 



da^ 



+ 6^" 



/h(") =pV|p|(a), (5.9) 



H,ff%{a, (/.) = I,. ;^|2 ^p(«> '^) • (5-10) 



with the term e^° acting as a potential barrier. At large negative a the solutions of this 
equation are oscillatory, whereas at large positive a the solutions are growing and decaying. If 
one were solving a scattering problem in Schrodinger quantum mechanics, one would demand 
that the solution decay at large positive a, in which case one would find equal admixtures of 
the oscillating solutions, e~*l^l° and e"*"*'^'", as a — > — oo. In our case, however, we demand that 
the solution behave as e"*'^'" as a — > — oo, in which case the growing solution is present (and 
dominates) as a ^ +oo. 

The dynamical behavior of <I> in the classically forbidden region e° ^ \p\ is most easily 
examined in the L^-version of the theory described in Sec. y. By Eq. (|2.1S| ), the effective 
Hamiltonian which governs the a-evolution of the L^-wavefunction, ^p(a, 0), corresponding to 
^p{a,(j)) is 

K 

I/h(«)P 

Since the function f\p\{a) grows rapidly for large a, the effective Hamiltonian approaches zero 
rapidly. This implies that, as claimed above, the dynamical evolution of (j) is "frozen" in the 
classically forbidden region, i.e., the probability distribution for (/> does not change with "time", 
a, when e^" ^ p^. 

In the region where e^" <C p'^, the wavefunctions <I>p(a,(/)) satisfy the free Klein-Gordon 
equation to an excellent approximation. Hence, there is no difficulty in constructing wavepackets 
which closely follow the classical trajectory in this region. On the other hand, as we have just 
seen, highly nonclassical behavior occurs for e^" ^ p^. It is of interest to examine the behavior 
of wavepackets in the classically allowed region, e^" < p^, when e^" ~ p'^ in order to determine 
more precisely where the semiclassical behavior breaks down. To investigate this, we consider 
the WKB approximation ||2^ and analyze the conditions under which it is valid and gives 
evolution corresponding closely to the classical dynamics. Since the equation we are solving 
is equivalent to the Schrodinger equation (^^), we see immediately that the WKB solution to 
order H which is positive frequency in the asymptotic past is 



•WKB/ - I 1*^1 * 



^H (") = vm^:^""p ~hi ^^"'^)^" 




(5.11) 
y Kyvt, P) I II J 

where 

' (5.12) 
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The corresponding L^-wavefunction is 



^ 



WKB 
P 



{a,(p) 



exp 



i / k{a,p)da + ip(f> 



(5.13) 



The criterion for the vahdity of the WKB approximation is given by |3C] 



2k{a,pY 



d_ 
da 



k{a,p) 



„Aa 



[p2 



4q,i3/2 



<1. 



(5.14) 



This inequahty is satisfied if p^ — e'^" ^ p^'^, which, in turn, is satisfied if p ^ 1 and p^'^ — e'^ ^ 
p~^'^ . Restoring the Planck length, /p, we findp ~ [Im/IpY , where Im is the maximum classical 
radius of the universe. Hence, the WKB approximation is valid provided only that Im ^ Ip and 
Im — l{oi) ^ [Ip/ImY'^Ipj where l{a) = e^lp. Thus, as long as the universe expands classically 
to a radius much larger than the Planck length, the WKB approximation is valid essentially up 
to the classical turning point. 

When the WKB approximation holds, the effective Hamiltonian is given by Eqs. ( |5.10| ) and 
( 5. lip . On the other hand, the classical deparametrized Hamiltonian (see Sec. y) with respect 
to the time variable a is simply Hd = k{a,p). Hence, when the WKB approximation is valid, 
we have Hgff ss Hd- Consequently, there should be no difficulty in constructing wavepackets 
which follow closely the classical trajectories pl| . 

The above conclusions can be verified numerically. Fig. 1 shows the quantity dcp/da = 
dH(.ff/dp as a function of a for p = 4, 16, 64 and 256. It is compared with the "true" classical 
trajectory, cp^i, given by Eq. ( [5.21 ). One can clearly see that the quantity dcp/da is approximated 
by d(l)ci/da better and better as p increases. Note also that the freezing of dynamics occurs 
more and more sharply as \p\ increases. Fig. 2 shows the wavepacket which, for a -^ — oo, takes 
the form 



$ 



p,Ap 



/2ApV/2 



\V^J 



exp <^ -2{Apy 



+ -ln2po 



tpo[a 



(5.15) 



with pq = 64 and Ap = 10. The shaded area represents the interval in (f) for each a which 
contains 90% of the squared L^-normalized wavefunction. This is compared with the classical 
solution a = (l/2)[lnpo ~ lncosh(2</>)]. This wavepacket indeed follows the classical trajectory 
quite closely up to the classical turning point at {(p,a) = (0, (1/2) Inpo)- Then it is seen to 
freeze near (j) = for a > (1/2) InpQ. 

Finally, we comment briefly on some differences between the approach taken here and those 
taken by Hartle and Hawking [ p2| , 33 1 and Vilenkin |34, |3^, |3^]. Our approach defines an 
entire Hilbert space of states, TC, and no state vector, ^ £ TC, is, in any sense "preferred". 
On the other hand, both Hartle and Hawking and Vilenkin seek to single out a particular 
wavefunction (without attempting to define a Hilbert space of states) via the imposition of 
boundary conditions. Nevertheless, one may inquire as to whether the Hartle-Hawking or 
Vilenkin wavefunctions lie in our Hilbert space, 7i. Even in the context of minisuperspace 
models, it is not clear how to implement, in a mathematically precise manner, the general 
principles which have been proposed to determine the Hartle-Hawking and Vilenkin boundary 
conditions. However, in the context of the simple model considered in this section, it seems likely 
that the Hartle-Hawking wavefunction is the solution to Eq. ( |5.3D for which $ is independent 
of (j) and $(a) — > const. / for a ^ — oo, i.e., ^HH{ct) = /o(e^"/2)- Similarly, the Vilenkin 
wavefunction should be the solution to Eq. ( |5.3| ) for which $ is independent of (p and <l>(a) 
decays for a — > +oo, i.e., $i/(a) = Ero(e^"/2). In any case, since the field transformation 
(p -^ (/i+const. is a symmetry of our model, it is clear that both the Hartle-Hawking and Vilenkin 
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wavefunctions (if they exist) should be independent of (p. On the other hand, any wavefunction 
in our Hilbert space approaches a Klein-Gordon normalizable free-particle wavefunction for 
a -^ — oo and the average of vr^ is necessarily nonzero. Hence, neither the Hartle-Hawking nor 
Vilenkin wavefunctions lie in our Hilbert space. 

6 The Taub model 



In this section we study the Taub model [37|, which can be obtained from the Bianchi IX model 
by letting /3_ = 7r_ = 0. From ([4.9| ) we see that the classical super-Hamiltonian for the Taub 
model is 

^T = -vr^ + 7r2^+e^"y(/3+)=0, (6.1) 

where 

1 _««, 4 



V{P, 



-e 
3 



-8/34 



-e 
3 



-2/34 



(6.2) 



As in the Robertson- Walker model of the previous section, all of the classical solutions start at 
an initial singularity at a ^ — oo, expand to a maximum size, and then recollapse. 
The super-Hamiltonian (|6.ip is simplified by defining 



r = (2a-/3+)/\/3, (6.3) 

e = (2/?+ - a)/V3 (6.4) 

so that the new coordinates, (t, ^), on minisuperspace differ from (a,/3+) by a Lorentz boost. 
In terms of the new coordinates, the Wheeler-DeWitt equation takes the form 

52 4 



Ht^ 



9r2 3 Oe 3 



$ = 0. 



(6.5) 



This equation can be solved exactly ||3^ by separation of variables. The equations for S^ and r 
both take the form 

"^ ■ "- F{x) = -oj''F{x) . (6.6) 



dx"^ 



Ae 



Bx 



The general solution to this equation is 



F{x) 



'2iuj/B 



'2^/A s. 



B 



-6 2- 



(6.7) 



where Ziy{z) is any modified Bessel function. The relevant solutions for TC are those of the form 

1 



^Ur,0 



Ur)SUO^ 



(6.8) 



where /(^(t) -^ e ^^'^ {lo > 0) for r ^ — oo and where Si^{^) has nonsingular ^-dependence. 
These conditions determine /^ and S^^ to be 



Ur) 



■)—iu)/\/li 



T{l-iu^/V?>)I_,, 



/V3U 



-e 



Vs? 



and 



s„(0 
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wsinh 



-KUJ 



/(2^/3) 



^/3^2 






(6.9) 



(6.10) 
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Here we have normalized S^j by requiring [p9| 



+00 



d(SU0Su.'i0 = si^-^') 



(6.11) 



and we have chosen the phase of Suj{C) to make it real. With this normalization, we have 

In this model, it is easiest to investigate dynamical evolution by choosing the "time slices" 
in minisuperspace to be the hypersurfaces of constant r rather than a. Again, we have A^ = 1 
and we choose iV" = 0. The L^-wavefunctions corresponding to ^cu{t,^) are given by 



^Ur,0 ^ Br<^Ur,0 = l¥rlsUO ■ 



and, by Eq. (2.19), the effective Hamiltonian is 



Heff^^{T,0 



LO 



\fUr)\ 



\Ur)\ 



r^.(r,0- 



(6.12) 



(6.13) 



Since fuiT) blows up for r -^ +cxd, we find that - as in the case of the scalar field dynamics in 
the Robertson- Walker model - Hgjj rapidly goes to zero, and the dynamics of (, with respect to 
r-time becomes "frozen" in the classically forbidden region, ^e 



4^2V3t 



>UJ' 



The WKB approximation can be considered for the r-dependence of the wavefunction with 
fixed CO. As in the Robertson- Walker model of the previous section, the criterion for the validity 
of the WKB approximation is satisfied in the classically allowed region, |e^^'^'^ < oj^, until very 
close to the classical turning point, provided only that the universe expands to a radius much 
larger than the Planck length. Thus, the (un-normalizable) basis functions for L^ states can be 
approximated by 



^(^(t-,0 ~ exp -i / dTK{T,uj) SujiO 



where 



k[t,u>) 



UJ 



4 2v^. 

3 



1/2 



(6.14) 



(6.15) 



Similar arguments apply to the function S^{^). In this case, the WKB approximation breaks 
down near the "spatial" classical turning point, i.e., for e~^^^^ ~ 3ca;^. However, as long as 
(jj^ ^ 1 , the region where the WKB approximation fails is small. Thus, except in a small 
region near the (spatial) classical turning point, the function Suj{^) can be approximated by a 
linear combination of the functions 



■^WKB 



(0 



UJ 



2TTk{^,L0) 



exp 



±i 



dmc,uj) 



(6.16) 



where 



k{^,u;) 



-4v^5 



1/3 



(6.17) 



Combining this result and that for the r-dependent part [Eq. ( |6.14| )], we conclude that a 
wavepacket can be made to follow the classical trajectory (except very near the point where 
the wavepacket bounces off the "wall" in ^) up to the classical turning point from expansion to 
contraction. 
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Again, this conclusion can be verified numerically. Fig. 3 shows a wavepacket with the 
central value of (D = w/yS = 32 and Au) = 5 which is arranged in such a way that it bounces 
off the wall before reaching the maximum expansion. It clearly shows that the wavepacket can 
follow the classical trajectory up to the classical turning point in r and then freeze there. 

The same analysis can be repeated for the Taub model with a homogeneous scalar field. 
The super-Hamiltonian of this model is given by 



Hts 



-7r„ + vr^^ + vr^ + e 



4a 



V{P, 



0. 



(6.18) 



One can readily solve the corresponding Wheeler-DeWitt equation by again introducing the 
variables r and ^, Eqs. ( |6.3| ) and (|67^). Since vr^ is conserved, we obtain the basis solutions 



^LopiT,i,i 



^/2t](ujj)J 



fri{uj,p){T)Sui{0 



2^ 



where ri{tj,p) = y/uj"^ + p^- The corresponding L basis functions, ^tjp, are given by 



^u>p{T,i,i 



and the effective Hamiltonian, -ffe/f ; is 



fr){uj,p) {t) 



fri{u,p){r) 



SUO 



Jp<t> 



(6.19) 



(6.20) 



-H'eff*a;p(-r,C,' 



r]{uj,p) 



fr,{LU,p) {t) 



r*L^p(T,^, 



(6.21) 



One advantage of this model is that there is more than one dynamical variable present, so 
that one could regard one of the variables - say, 4>, since classically it evolves monotonically - 
to be a "physical clock", against which the evolution of the other variable, ^, is to be measured. 
However, since fri(uj.p){T) — > oo as r ^ oo, it is clear that all of the dynamical variables, 
including the "physical clock", will "freeze" in the classically forbidden region. Thus, we see 
that the "stoppage of motion" predicted in this model actually would be physically unobservable; 
rather, this freezing would correspond, physically, to measurable time stopping at the maximum 
expansion of the universe. 

It would be of interest to investigate whether the above "freezing" phenomenon persists for 
other choices of time slicing - particularly the a-slicing - and, if so, which variables "freeze" . Our 
calculations above for the r-slicing were enormously simplified by the presence of an orthogonal 
basis of solutions in a separated form. Since the r-derivatives of these basis solutions are 
proportional to themselves, Ct is diagonal in this basis. Thus, one can easily obtain Cj, and, 
thereby, calculate Hgjf. On the other hand, for the a-slicing, although it is not difficult to 
calculate Ct directly from the definition ( |2.9p , we do not thereby obtain simple expression for 
CJ. For this reason, we have been unable to calculate H^ff and study the issue of "freezing" in 
the Q-slicing. 
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Appendix A: An alternative quantization of class A Bianchi mod- 
els 

In Sec. ^we eliminated the momentum constraint of the class A Bianchi models by noting that, 
classically, we can choose ha/s and vr"^ to be diagonal. We then wrote down the Wheeler-DeWitt 
equation for the diagonal Bianchi metrics and found it to have the form of a Klein-Gordon 
equation in a 3-dimensional flat spacetime with a time-dependent potential. However, it is 
more natural not to impose the diagonal form of the metric at the outset, and to proceed by 
treating the Hamilitonian and momentum constraints on an equal footing. We employ this 
latter approach in this Appendix. For simplicity, we specialize to the vacuum models and will 
focus upon the Bianchi I and Bianchi IX cases. Many of the results here, including the algebra 



of the momentum constraints, have been obtained by Ashtekar and Samuel |40|. 

We consider, first, the form of the Wheeler-DeWitt equation for the Bianchi models when the 
restriction to diagonal metrics is not imposed. For each Bianchi Lie group, G, minisuperspace is 
comprised of the 6-dimensional manifold of left invariant metrics, hat, on G. The Hamiltonian 
constraint is 

^ = G',fe,dvr"V'^ + F = 0, (Al) 

where the inverse supermetric is 

Gabcd = -^h^ {hachbd + had^bc — habhcd) (A2) 

and the determinant, h, of hab is defined with respect to a fixed, left invariant volume element 
on G. It should be noted that the two lower index pairs {ab) and {cd) correspond, respectively, 
to the upper indices A and B in Eq. ( |2.1| ) of Sec. |2[ Thus, the supermetric. Gab, is 



f-iabcd 



-h^l"^ {h'^^h^'^ + h'"^h^^ - 2/i"''/i^'^) . (A3) 



The geometry of Gab was studied in the classic work of DeWitt [^] . The signature of Gab 
is Lorentzian ( — h + + ++), and the transformations hob — > e'^°'hab define timelike, hypersurface 
orthogonal, conformal isometries of Gab- We may identify each of these orthogonal hypersur- 
faces with the 5-dimensional manifold, TWc of left invariant conformal metrics hab = h~^'^hab 
on G. The supermetric then takes the form. 

Gab = e=^" (-24VAaVBQ + Hab) • (A4) 

The Riemannian metric, Hab, on A^c" is invariant under the SL{3, IR) transformations 

Kb^Na'N.'^h.d, (A5) 

where N^'^ satisfies detA^ = 1. If we introduce arbitrary local coordinates ^ (A = 1 — 5) on 
A^c, then the coordinate components of Hab take the explicit form 



-^ab ^hbc Tcddhda 



HA^ = h^'l^h'''^, (A6) 



where /i"* is the inverse of hab- The Ricci tensor and the scalar curvature of Hab are 

Rab = -^Hab, (A7) 

R = —^ (A8) 
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(see Ref. pl[] , Appendix A, for details; we have corrected the corresponding equations of that 
reference by a factor of 1/2.) Thus, in contrast to the diagonal case, the geometry of mini- 
superspace is now curved. Thus, the Wheeler-DeWitt equation ( ^^ ) corresponding to (^ 
now takes the form of a curved space wave equation 



da^ 



IADaD"^ - 90^ + 241/ 



$ = , (A9) 



where Da is the derivative operator associated with Hab- 

For the Bianchi I model, no momentum constraints are present so the full quantum con- 
straints will be imposed on $ by satisfying Eq. (A9) with V = Q. Since Eq. ([A 9]) is a curved 



space equation, it is easily seen that the theory obtained in this manner is not, in any sense, 
equivalent to the theory obtained by restricting to diagonal metrics as done in Sec. |^. (The 
reasons for this inequivalence were explained in that section.) 

One further feature of Eq. ( [A9| ) in the Bianchi I model is worth noting. Since minisuperspace 
is six dimensional, the natural Q , conformally invariant version of the Wheeler-DeWitt equation 
is obtained by choosing ^ = 1/5. However, the lower bound of the spectrum of —DaD is 



1/2 (see, e.g., Ref. [42|, p. 49). Thus, the Wheeler-DeWitt equation with this choice of ^ 
is "tachyonic", and there exist spatially well behaved solutions which grow exponentially as 
a -^ iboo. In particular, the conformally invariant choice of Wheeler-Dewitt equation does 
not allow the division of the space of solutions into asymptotic positive and negative frequency 
subspaces as a — > — oo, i.e., our prescription for the construction of TC breaks down. However, 
this difficulty can be avoided by making a different choice of ^. 
The classical momentum constraints 

- 2Da1T'''' = (AlO) 

are nontrivial in all the other Bianchi models. (Here Da denotes the derivative operator on G 
associated with the left invariant metric hat-) In order to motivate a natural operator version 
of these constraints, we rewrite them as follows. First, if s" is a left invariant vector field on 
G, then by a direct computation using the formula for the connection in terms of the structure 
constant field G"^^ (see, e.g., Ref. [||], p. 314), we obtain Z^^s" = G^^s" = Aas"" = 0. Hence, 
we find that the momentum constraint 



can be rewritten as 



Ps = -2s'hcDa7r'''' (All) 

Ps = -2Da{s'hcTT''') + 2{DaSb)7r'''' 

= {£shabK\ (A12) 

since the first term on the right side of the first line of this equation is the divergence of a left 
invariant vector field and, thus, vanishes. 

When we promote Pg to an operator, it is natural to choose the factor ordering ( [A12| ), 
since then the constraint Ps^ = implies that the wavefunction $ is invariant under the 
transformation 6hab = Cghat- (As is well known, an analogous choice and interpretation of the 



momentum constraint operators can be made in full quantum gravity [41|, the only difference 
being that neither hat nor s" are restricted to be left invariant.) Note that with this choice of 
P,, we have 



[hab,Ps] = iCshabi (A13) 

= iCsT^''\ (A14) 



^"Kp. 
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and 



[Ps,Pt]=iPlsA 



(A15) 



Let us study the imposition of the momentum constraints for the case of the Bianchi IX 
model in some detail. The Lie group, G, of the Bianchi IX model is SU{2), and the inverse, 
{M~^)ab, of the left invariant tensor field M"-^ defined by Eq. ( |4.3| ) provides a positive definite 
metric on the 3-dimensional vector space, W, of left invariant vector fields. The exponenti- 
ated version of the momentum constraints P,^ = imply that the allowed wavefunctions are 
invariant under the transformation 



h 



ab 



OaOb hcd : 



where O^ is an arbitrary orthogonal linear map [with respect to {M^^)ab\ on W with positive 
determinant. 

It is convenient to introduce the following parametrization of minisuperspace. Since hab 
defines a positive self-adjoint map on W in the inner product (Af"^)^^, we may characterize 
each hab (i-e., point in minisuperspace) by its eigenvalues, Ai > A2 > A3 > 0, and the rotation, 
O, which carries its triad of eigenvectors into some fixed, orthonormal basis. Note, however, 
that this parametrization becomes singular when degeneracy occurs in the eigenvalues of hab- 

The Hilbert space of states, W, is constructed from functions <I>(Ai, A2, A3, 0) on mini- 
superspace which satisfy the momentum and Hamiltonian constraints. The momentum con- 
straints are easily imposed, since they simply require <I> to be independent of O. The form of 
the Wheeler-DeWitt equation in these coordinates can be obtained as follows. On the manifold, 
A^C", of conformal metrics (where A1A2A3 = 1, so only two of the Aj's are independent), we have 
from Eq. ([A^) 



^AS 



Tr 



h-'^h-' 9/^ 



-Tr 



' dT dT-^\ 



2Tr 






o 



-'''-TO-^S.T- 



dC 



dc} 



(AI6) 



where Tab is the diagonal matrix with eigenvalues Ai, A2, A3. By writing 

O-^dO = 






-^3 


X2 


^3 





-Xi 


-X2 


Xi 






(A17) 



we obtain |42| 



+Ai(A2 - Xsfxf + A2(A3 - Xifxl + A3(Ai - X2fxl 



(A18) 



Although the momentum constraints require ^ not to depend upon O, the portion of the 
supermetric involving (Xi, X2, X3) makes a A-dependent contribution to the volume element on 
superspace, and this, in turn, nontrivially affects the form of Wheeler-DeWitt equation. Taking 
into account the condition A1A2A3 = 1 on Mci we see that this portion of the metric contributes 
to \/H the additional factor 



C 



n(A. 



A,) 



(A19) 
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Writing Ai = e2(/5++v^/3-), A2 = e2(/3+-v^/3-), and A3 



-4/34 



, we obtain 



C{I3±) -- 
The Wheeler-DeWitt equation then becomes 



sinh2\/3/3_ sinh(3/3+ - \/3/3_) sinh(3/3+ + \/3/3_ 



' E|;^^(/?±)|;r + 24F-9oe 



aa2 C(/3±) ,tl 9A 



aA 



$ = 0. 



(A20) 



(A21) 



The variables 13± are restricted by the condition Ai > A2 > A3 > 0, but we may drop this 
restriction by requiring, instead, the invariance of ^ under the transformations of I3± corre- 
sponding to permutations of Aj (i = 1, 2, 3). The kinetic term in the Wheeler-DeWitt equation 
can be converted into flat spacetime wave operator by rescaling <1> as <1> = C(/3±)~^"$. Then, 
we obtain 

+ 241/ - 90e 



da^ dpi dpi 



with 



V = V + 



24C7(/3±)V2 



^2 



+ 



92 



dpi dP: 



$ = 



C(/3±)V2 



(A22) 



(A23) 



This differs substantially from the form of the Wheeler-DeWitt equation derived in Sec. ^. 
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Figure captions 

Figure 1: d(l)/da = dHejj/dp for p = 4, 16, 64, 256. The dotted line is the corresponding classical 
value [1 - e4°/p2]-i/2. 

Figure 2: The region swept by the interval of (j) (for given a) containing 90% of the squared 
L^ wavefunction for a Gaussian wavepacket with (p) = 64 and Ap = [{p"^) — (p)^]^'^ = 10 in 
the Robertson- Walker model with a massless scalar field. The solid line is the corresponding 
classical solution. 

Figure 3: The region swept by the interval of ^ containing 90% of the squred L^ wavefunction 
for a Gaussian wavepacket with {to) = 32 and Atu = 5 in the Taub model, where we have defined 
(D = 3~^''^L0, ^ = 3^'^.^, f = 3^' ^r. The solid line is the corresponding classical solution. 
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